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Abstract. Within the framework of a flat cosmological model a propagation of an instanta- 
neous burst of nonpolarized isotropic radiation is considered from the moment of its beginning 
at some initial redshift zq to the moment of its registration now (at z = 0). Thomson (Rayleigh) 
scattering by free electrons is considered as the only source of opacity. Spatial distributions of 
the mean (over directions) radiation intensity are calculated as well as angular distributions of 
radiation intensity and polarization at some different distanses from the center of the burst. 
It is shown that for redshifts zq large enough (zq > 1400) the profile of the mean intensity 
normalized to the total number of photons emitted during the burst weakly depends on initial 
conditions (say the moment z of the burst, the width and shape of initial radiation distribu- 
tion in space). As regards angular distributions of intensity and polarization they turn to be 
rather narrow (3 - 10 arcmin) while polarization can reach 70%. On the average an expected 
\ polarization can be about 15%. 

Key words: cosmology, early Universe, cosmological recombination, radiative transfer, Thomson 
scattering. 
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^2 ' INTRODUCTION 

Investigations of cosmic microwave background (CMB) both theoretical and experimental 
ones continue extensively all over the world. Achivements on this way are enormous. Infor- 
mation about fundamental parameters of the Universe is obtained. Contribution of completely 
new kind of matter inaccesible for investigation in ground based laboratories is discovered and 
measured. The main conclusions of the theory of origin of large scale structure of visible mat- 
ter distribution in the Universe are confirmed. However the prospects of future investigations 
seem to be still more grandiose ones. In particular completely new discoveries are possible in 
the course of thorough studying of small scale CMB temperature fluctuations. Our paper is 
devoted to one of the mechanisms of such fluctuations. 

In standard cosmological model fluctuations of matter density are possible on different 
scales - from very large to comparatively small ones. An interaction of these fluctuations with 
CMB during Universe expansion and cooling determines observed pattern of CMB temperature 
variations (AT/T) at the moment of hydrogen recombination at a redshift z fh 1000. Scales 
smaller then someones "damp" i.e. give very small values of AT/T. This takes place for different 
reasons, namely through decrease of peculiar velocity of matter density fluctuation (cloud) due 
to radiation viscosity which leads to decrease of a Doppler shift of photons frequency; through 
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decrease of cloud size relative to horizon and decrease of contribtution of Sacks- Wolfe effect; 
through decrease of cloud optical thickness for Thomson scattering due to electron density 
decrease. However decrease of amplitude of the own matter density fluctuation takes place in 
this case for sonic waves only! At the same time for some class of fluctuations their contrast 
relative to surrounding matter increases. All these results are true for standard small density 
fluctuations. But we consider an evolution of AT/T for nonstandard models of origin of energy 
release on small scales. In particular this can be primary black holes, domains of different 
unstable forms of matter etc. at the epoch before hydrogen recombination. In such an aspect 
this problem was formulated in the paper by Dubrovich (2003). In the same paper an evolution 
of spatial distribution of "superequilibrium" radiation intensity was described qualitatively. 
Conclusion was obtained that an angular size of AT/T fluctuation observed now did not depend 
on a redshift of radiation source. 

The aim of the present paper is a correct calculation of this process within the framework 
of standard theory of radiative transfer. As a model we use a source of sufficiently small size 
radiating during some small interval of time (instantaneously in particular). We assume that 
this radiation does not have any influence on the matter parameters (e.g. electron temperature 
and number density) during the whole process of scattering. Within the framework of the 
present paper we do not discuss an absolute quantity of scattered radiation intensity as well 
as its spectrum. Separate paper will be devoted to this question. We stress once more that 
we consider here completely different type of temperature inhomogeneity as compared with the 
standard scenario. Our source of inhomogeneity is not a small fluctuation of matter density 
described in the course of its evolution by equations of hydrodynamics and thermodynamics. 
Therefore conclusions of the standard theory of small scale fluctuations damping (Silk damping) 
are not applicable to it. 



We consider nonstationary radiation transfer in homogeneous expanding Universe starting 
from spherically-symmetrical initial distribution of radiation intensity. For a flat model of the 
Universe the corresponding equation of radiation transfer for photon occupation number n has 
the form (Nagirner, Kirusheva, 2005): 



where rj is a conformal time (drj = cdt/a(r]), c is the speed of light, t is a time, a = a(rj) is 
a scale factor), p is a spatial coordinate, H = H(n) is a Hubble factor, n = (c 2 /2hv 3 )(I, Q) T , 
I = I(p, ii,r],b>) and Q = Q(p, p, r], v) are Stokes parameters for radiation at a frequency v 
propagating at an angle $ = arccos p to radial direction, dimensionless vector source function 
is s = (si, sq) t = (c 2 /2/iz/ 3 )(S7, Sq) t = (c 2 /2hu 3 )e/k and e = e(p, p, 77, z/) and k = k(p, p, r], v) 
are emission and absorption coefficients accordingly. It should be noted that k = k{rf) for 
Thomson scattering in a homogeneous medium. Here and further T means transpose. 

Below we assume scale factor a to be equal 1 at the present epoch (z = 0) so that 
a(z) = 1/(1 + z) where z is a redshift. Therefore spatial coordinate p in eq. (pQ) is a distance 
from a center of symmetry measured in the present epoch (at z = 0). At an arbitrary z the cor- 
responding distance is r = pa(z) — p/(l + z). Further since Thomson scattering is neutral one 
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than radiation frequency v changes due to cosmological redshift only: v = v^jaiz) = u Q {l + z), 
where u is a frequency at the present epoch. Therefore we may use u instead of v as a 
frequency variable and we redefine the source function (Nagirner, Kirusheva, 2010) as follows: 

n(p, fi, r], v) = n*(p, p, rj, u ) = n ( p, p, 77). (2) 

Then frequency variable disappears and eq. ([1]) takes the form 

dfi dn 1 — p^ dfi 

drj + ^~dp + p dp = ~ k ^ a ^ /i ' ^ ~ /i ' ^ ' ^ 

In the case of monochromatic scattering and in the presence of axial symmetry of radiation 
field the vector source function in the righthand side of this equation is written in the form 

s(p,V,v) = (V 2 ) / d(iR(ji,fi)fi{p,(i,ri), (4) 



where phase matrix R(n, //) describes radiation redistribution over directions and polarizations 
in a single scattering. For Thomson (Rayleigh) scattering phase matrix can be factorized as 
follows (see e.g. Ivanov, 1995): 

i?(//,//) = i(/i)i T a/), (5) 

where 



In this case eq. (jl]) takes the form 

s(p,V,v) = A(ji)s*{p,ri), (7) 

where 

s*(p, V ) = (1/2) A T (p')n(p,p', V )dp'. (8) 

Eq. ([2]) looks like an ordinary nonstationary equation of monochromatic radiation transfer 
in a medium with the volume absorption coefficient 

a( v ) = k(rj)a(ri), (9) 

which depends on time only. Therefore 77/c can be considered as an ordinary time and p as 
an ordinary radial distance in the problems of nonstationary radiative transfer in spherically- 
symmetrical extended media. Thus we can use the method suggested earlier by Ivanov (1995) 
for solving problems of stationary transfer in plane-parallel media with generalized Rayleigh 
scattering. The method consists in reducing the problem to getting and solving an integral 
equation for so called reduced source function defined by eq. (jSJ) which does not depend on 
angle variable. After that radiation field with its angular structure is defined from the formal 
(i.e. with the known source function) solution of initial eq. 

To obtain above mentioned formal solution we consider propagation of radiation along a 
ray intersecting radial direction at an angle = arccos p, at a distance p from the center of 
symmetry. Let us introduce coordinate I measured along a ray from the point nearest to the 
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center of symmetry in direction of radiation propagation. The lefthand side of eq. ([3]) is nothing 
else than derivative along the ray direction so that this equation is rewritten as 

dn/dl = -a(l)[n(l)-s(l)}. (10) 

Integrating this equation we obtain 

n(l) = n (h)e' J 'o ail ' )dl ' + t a(l')s(l')e~ fi a ^ dl " dl' . (11) 

Jlo 

Here Iq is a coordinate of a point farthest from the observation point I but yet capable of making 
a contribution to radiation in this point at a given time r]. Each point /' on a ray is defined 
by radial distance p' and by an angle arccos p! at which the ray intersects radial direction. It 
follows from geometry of the problem that 

I = pp, I' = pp, l = p po, p\Jl- p 2 = p'y/l - p' 2 (12) 

and 

l-l = r), I -I' = rj-rj'. (13) 

According to the last of these equations one can turn to integration over time in eq. ( TTTj) 
because dl 1 = drj'. As a result the formal solution (TTTj) can be written as follows 

n(p,p, V ) = no(p ,p )e-fo"W'W + f a (rf)A(ji'W, i)e~ § aW ' W 'drf, (14) 

Jo 

taking into account eqs. (jTJ), f|T2j) and (TT~3]) . Here vector n (p, p) defines initial (at the moment 
rj — 0) distribution over distances and angles and 

po = \jp 2 ~ 2ppr] + r] 2 , popo = pp-r], (15) 

P = \jp 2 ~ ZpfMjl - V') + (V- V') 2 , p'rf = pp-V + v'- (!6) 
It follows from the formal solution f)14p that one can introduce dimensionless time 

u = [ a{rf)drf = ca e [ n e (t')dt' = co e [ — dz' (17) 
Jo Jo Jz (1 + z )H \z ) 

which has a sense of an optical distance (by Thomson scattering) between moments z and z . 
Here redshift Zq corresponds to the initial moment of time: u = rj = t = for z = Zq, n e is 
an electron number density, a e = 6.65 • 10~ 25 2 is the cross-section of Thomson scattering. For 
conformal time r] we have equation 

rt rzo 
7] = c dt'/a{t') = c / dz'/H(z'), (18) 

Jo Jz 

which can be used to relate u with rj by calculation both of them on the same grid of redshifts 
z. With the new time variable u the formal solution (Tl4"]) takes the form 

n(p,p,u) = n (p ,po)e~ u + / A(p')sjp' \u')e u ~ u du' . (19) 

Jo 
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Substitution of eq. (|T9|) into the righthand side of eq. (|8]) gives the following equation for 
s*(p,u) = s (p,u) + - / e u ~ u / A T (p)A(p')s4p f ,u')p , dp f , (20) 

2p J0 1] - rj J\p-r,+r,'\ 

where arguments of matrices A T and A are equal (according to eq. ffl6|) ) to 

a* = [p 2 - p 2 + (v- vTUiMv - v% p = [p 2 - P 2 -(v- v')W(v - v% (21) 

Further the primary source vector in eq. fl20l) is 

e~ u rp+v 

Mp> u )=o — / A (p)no(Po,Po)Podp , (22) 

ipj] J\p-r]\ 

where 

P=(p 2 -P 2 + V 2 )/2PV, Po = (p 2 ~pl-V 2 )/2poV- (23) 
according to eq. f|T5l) . 

To deduce the main integral equation fTSU]) we passed from integration over p to integration 
over p' in the integral term and to po i n the free term using in the first case the first of eqs. 
(I2~T|) which gives 

dp = -p'dp'/p(r] - r/), (24) 
and in the second case we use the first of eqs. fT2"3l which gives 

dp = -podpo/pr]. (25) 

In scalar case assuming isotropic scattering the main integral equation (TSUI) has a more 
simple form 

1 r u 1 dv! rp+v( u )-v( u ') 
s(p,u) = s (p,u) + — e u ~ u —- — r / s(p,u')p'dp, (26) 

2p Jo 7?(W) - r]{u') J\p- v (u)+rj(u')\ 

where 

s (p,u) = - — / n Q (po)podp , (27) 



2pr] J\p- V \ 

and from eq. (T191) we obtain for the formal solution: 

n(p, p, u) = n (p )e~ u + T s(p', u')e u '- u du', (28) 

Jo 

where 

Po = \Jp 2 -2ppr] + r] 2 , p = \J p 2 - 2pp(r] - rf) + (77 - 77') 2 - (29) 

It should be stressed that in this case sourse function coinsides with the mean (over angle 
variable) radiation intensity: 

s(p,u) = j(p, u) = (1/2) f n(p,p,u)dp (30) 
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Here and below the term "intensity" means dimensionless intensity i.e. photon occupation 
number. 

As an initial condition we assume that at the moment t = (rj = 0) corresponding to some 
redshift z$ radiation is nonpolarized and isotropic and has spherically-symmetrical distribution: 

n(p,p,0) = K(p),0) T , (31) 

where no(p) is a given function which we take in the form 

n (p) =7r- 3 / 2 p; 3 exp[-(p/p*) 2 ] ^5(p)/(4np 2 ) for p* -> 0, (32) 

where p* is a parameter wich can be taken sufficiently small to model point source. Obviously 
no (p) satisfys to normalization 

poo 

4tt / n (p)p 2 dp = 1. (33) 
Jo 

We solve the problem of propagation of instantaneous burst of radiation in a scattering 
expanding and recombining Universe. Since albedo of Thomson scattering is equal to 1 the full 
number of photons emitted in the burst must be conserved. So that an equality 

roc /"OO 

4tt / j(p, u)p 2 dp = 4tt / n (p)p 2 dp = 1. (34) 
Jo Jo 

must be fulfilled at any moment of time. One can make certain that this equality follows indeed 
both from vector equation (120]) and scalar eqution ([26]) . 

METHOD OF SOLUTION AND MAIN RESULTS 



We obtained numerical solutions of vector equation ( 1201) and scalar equation ( |26|) by means 
of their discretization on given grids over dimensionless time u and over distance p measured 
in Mpc. It should be noted that the scale of distances p corresponds to the present epoch 
according to our normalization of scale factor: a = 1 at z = 0. As a main time grid we use 
a uniform grid over redshift z with the step Az = 10 and grids over u and rj are calculated 
then using eqs. lITTj) and (|T8l) . The function s*(p', v!) in the righthand side of integral equation 
(1201) was approximated as a function of p' for fixed u' by cubic spline on a uniform grid with 
the step Ap = 1 Mpc and integral over p' was calculated analytically. Next the whole of 
integrand (except for exponential factor) in integral over u' was approximated (as a function 
of v!) by quadratic spline and integral over u' was calculated analytically as well. At the last 
time point (for u' = u, rj = rj) the whole integrand in integral over u' including the multiplyer 
l/2p turns to be equal to Bs*(p, u) where B =diag(l,7/10) is a diagonal matrix. Carrying this 
term of quadrature sum from the righthand side of equation to the lefthand one we get in the 
end recurrence relation which allows to express the current solution through solutions at the 
preceeding moments of time. To control the process of solution we check up conservation of 
the full number of photons. It was fulfilled with a relative error no more than 1CT 7 for scalar 
equation and no more than 10 -5 for vector equation. 

After the source function was found then angular distributions of radiation intensity and 
polarization were calculated numerically from the formal solution for different distances from the 
center of the burst and at different moments (z) including the present one (z — 0). Calculations 



6 



were fulfilled for several values of initial moments of time in the range of redshifts z from 1100 
up to 3000. 

The width of initial intensity distribution as a function of p (see eq. (152]) ) was taken to be 
p* = 1.5 Mpc in the scale of distances at z — 0. But in the scale of distances corresponding to 
the moment of the burst (at some z = Zq) the width of initial distribution becomes (for zq ^> 1) 
much smaller: r* = a(zo)p* = P*/(l + Zq)- 

As regards to another parameters appearing in the problem they enter the Hubble factor in 
particular: 

h(z) = Ho^/n A + (i - + zf + n M (i + zf + n rel (i + z)\ (35) 

where H = 2.4306 • 10~ 18 h s -1 , ho is the Hubble constant in the units of 75 km/ (s-Mpc); Qm, 
Qa and Q re i are ratios of densities of matter, dark energy and relativistic particles (radiation, 
massless neutrino) to the crytical density p c = 3Hq/(8ttG) at the present epoch; Q = Qm + 
Q\ + Q re i, &rei = Pj?(l + fn)/Pc Pr = urTq/c 2 is a mass density of radiation at the present 
epoch (To is the mean temperature of CMBR), f n is a contribution of relativistic (massless) 
neutrino (usually f n = 0.68). For the flat model of the Universe we have Q = 1 and then 

Moreover number density of electrons enter the equations. It is measured usually in the 
units of the total number density riu of hydrogen atoms and ions: n e (z) = x e (z)nn(z), where 
x e (z) is so called recombination history of the Universe and 

n u (z) =n° u (l + z) 3 , < = 0.63144 • 10" 5 Xfi B ^o " 3 , ( 36 ) 

where is a ratio of barion density to critical density at the present epoch, X is a hydrogen 
abundance (by mass). Recombination history is calculated separately and enter as an input 
file. We calculated it using the programme "recfast.for" (Seager et al., 1999). 

We adopted the following values of parameters: Q — 1, Q\ — 0.7, Qb = 0.04, T = 2.728 K, 
hydrogen abundance X = 0.76, Q re i = 0.85 • 10 -4 , Hubble constant H = 70 km/(s-Mpc). 

Fig. 1 shows initial isotropic distribution of radiation. Results of calculations are shown in 
Figs. 2-7. Figs. 2-5 are obtained in scalar approximation and Figs. 6 and 7 - using exact 
description of scattering. Fig. 2 displays propagation of radiation wave generated by the burst 
at the epoch (z = 10 and 20) when the Universe was practically transparent for radiation (we 
do not consider here reionization connected with the birth of primary stars). Then radiation 
propagates freely and an interval between peaks in Fig. 2 is equal to the difference of conformal 
times of the burst for corresponding redshifts. 

In Fig. 3 there are shown distributions of the mean intensity (i.e. source function) at the 
present epoch but for different moments (zq) of the burst. Certainly all profiles are in the 
range of distances which do not exceed (strictly do not exceed if we neglect the width of initial 
spatial profile of the burst) conformal time of the burst (see Tabl. 1) since according to its 
definition (eq. (I18p) conformal time is equal to a distance passed by freely flying photon from 
the moment of its emission (at some Zq) to the moment of its registration (at z = 0) in a fixed 
scale of distances coinsiding by our choice with the scale in the present epoch. Also it is clear 
that distribution of diffuse radiation should be wider and its maximum should be nearer to the 
center of the burst as compared with straightly passed unscattered radiation which should have 
its maximum at a distance equal to conformal time of the burst. While maximum of diffuse 
radiation is situated according to our calculations at a distance p = 13620 Mpc irrespective of 
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the moment z of burst beginning. This distance coincides as it is expected with the conformal 
time for the last scattering surface at z ~ 1090. This is illustrated by two-component profile 
in Fig. 3 for zq = 1200 containing both diffuse and straghtly passed radiation. 

Furthemore it is seen in Fig. 3 that for sufficiently large zq > 2000 the profile of distribution 
does not depend practically on the time of the burst beginning. Moreover our calculations show 
that it does not depend (providing that the full number of photons conserves) on a characteristic 
size of the burst (in our case this is parameter p* in the initial distribution ( 1321) ) provided it 
is not too large (see below). So for sufficiently large zq the profiles of intensity distribution 
should be in a sence universal ones i.e. nondepending on initial conditions of the burst as it 
was indicated earlier by Dubrovich (2003). This is connected with the large optical thickness of 
the Universe by Thomson scattering at large redshifts (see Tabl. 1). So that photons emitted 
during the burst are trapped and a size of radiating region changes due to diffusion not too 
much till the moment of sufficient clearing of the Universe at z ~ 1200 owing to hydrogen 
recombination. 

It should be noted that in accordance with our solution normalization on the full number 
of photons emitted during the burst into space the calculated quantities have a dimension of 
inverse volume, namely 1 Mpc~ 3 . Therefore if at the moment of the burst its central intensity 
(at p = 0) is equal to n then our profiles must be multiplied by 7r 3 / 2 n (p*/lMpc) 3 according 
to the form of our initial profile (eq. f[3"2"j) ). But in general case taking into account weak 
dependence of solution on the initial conditions one can take as a transitional multiplier e.g. 
(47r/3)n(p /lMpc) 3 where p is a typical size of the burst and n is intensity averaged in a sphere 
of this size. 

In Fig. 3 one can see also that as redshift zq approachs to the beginning of hydrogen recombi- 
nation the changes of the profile become more and more appreciable and during recombination 
when clearing of the medium becomes noticable the straightly passed unscattered radiation of 
the burst appears: at z Q = 1200 it is seen as a narrow peak to the right of diffuse maximum 
and at z = 1100 it is already dominates (see Fig. 4). 

In Fig. 5 one can see angular distributions of radiation intensity at different distances from 
the center of the burst and at different moments (zq) of the burst (see Tabl. 2). Characteristic 
feature of these distributions is their small width (3 - 10 arcminutes) and the width is smaller 
on the leading front than on the rear one and it decreases with the growth of Zq. It should be 
stressed that owing to such a strong anisotropy the radiation intensity towards the center of 
the burst can be almost 10 7 times larger than angle averaged intensity (see Tabl. 2 and Fig. 
3). It is clear that degree of diffuse radiation anisotropy is defined by an angle 2d d at which 
radiating region is seen at the present epoch (z = 0) on the last scattering surface (z\ s « 1090, 
77i s ~ 13620 Mpc). If r d is the region radius at that moment then obviously /i^ = cos$d ~ r d /rji s 
for r d = PdjiX + Zis) <C 77is. For straghtly passed radiation the cosine of angle at which the 
burst (arised at z = z ) is seen at a distance p from the center amounts p r = cos^ r ~ r*/p for 
= p t /(l + zo) C where p* is the initial radius of the burst in the present scale of distances. 
For example for the burst at zq = 1100 we have p = 77 « 1.36 • 10 4 Mpc taking p* = 1.5 Mpc 
for the initial radius of the burst. So that we have p r ~ 10 -7 and $ r ~ 1.5 arcminutes. In this 
case as was already pointed above the main part of radiations comes to us without scattering 
on the way. As for diffuse radiation which dominates for z Q > 1200 the width of its angular 
distribution turns out to be larger. Thus for zq > 1600 the semiwidth of angular distribution of 
radiating region at z = is equal to 3 - 10 arcminutes according to our calculations so that its 
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u 


i], Mpc 


3000 


177.23 


13770 


2000 


67.63 


13720 


1600 


29.27 


13686 


1400 


12.74 


13664 


1200 


3.01 


13635 


1100 


1.07 


13618 



Table 1: Optical distances u and conformal times rj of the burst from its beginning at z to the 
present epoch z = 0. 



N 


r, Mpc 


z = 1600 


z = 2000 


z = 3000 


z = 2000 


1 
2 
3 
4 


13570 
13600 
13630 
13660 


1.02 • 10- 6 

6.91 • 10~ 6 

1.92 • 10~ 5 
6.97- 10~ 7 


1.02 • 10~ 6 
6.56 • 10~ 6 
1.57- 10~ 5 
7.89 • 10~ 7 


1.02 • 10- 6 
6.28 • 10~ 6 
1.38 • 10~ 5 
8.60 • 10~ 7 


1.27- 10~ 6 
7.39 • 10~ 6 
1.69 • 10~ 5 
9.00 • 10~ 7 



Table 2: Radiation intensity n(0) towards the center of the burst as a function of distance r 
from the burst center for different moments zq of the burst beginning. Last column corresponds 
to exact description of Thomson (Raylegh) scattering and preceeding three columns correspond 
to scalar approximation. 

radius at z — 1100 must be 2 - 7 times larger than for the burst at z = 1100 i.e. p d 3 — 10 
Mpc in the present scale of distances and r d = p<i/(l + z ) in the scale of distances at z — z . 
Therefore the characteristic initial radius of the burst at z > 1600 must be in any case smaller 
than this estimate lest it should influence vitally on the propeties of the burst radiation at the 
present epoch. 

Further, decrease of anisotropy when passing from the "base" of leading wave front over 
maximum to the "base" of rear front (see Fig. 3 and Tabl. 2) is explained as follows: photons 
observed at the largest distances come from the nearest part of radiating region in a small solid 
angle and an effective size of region giving a contribution to observed radiation grows with the 
distance decrease. 

Fig. 6 shows change of angular distributions when passing from approximate scalar descrip- 
tion of Thomson scattering to an exact description taking into account scattering anisotropy 
and polarization. Finally Fig. 7 shows that polarization of radiation can be rather large (up to 
70%) but on angular distances where radiation intensity becomes already much smaller than 
towards the center of the burst (see preceeding Figure). It should be noted also anticorrela- 
tion between polarization and anisotropy of radiation. Namely, when passing from the farthest 
(from the burst center) point of the wave profile to the nearest one an anisotropy decreases 
but polarization grows (cf. Figs. 6 and 7). The thing is that radiation from farthest points 
comes from the nearest to observer small part of emitting region and mainly it consists of pho- 
tons undergone their last scattering almost directly forward which changes polarization only 
slightly. But for the less distant (from the burst center) points an essential contribution comes 
from photons scattered under sufficiently large angles which leads to polarization growth. 
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CONCLUSIONS 



A source of very small size radiating energy in prerecombination epoch will be seen now as 
a some spot on the background of cosmic microwave radiation. Our calculations confirm initial 
conclusion made by Dubrovich (2003) about nondependence of angular size of this spot on the 
moment of the source burst on condition that it is situated at a distance of optimal visibility. For 
estimates of radiation intensity it is important that duration of the source burst also has a very 
little influence on the size of the spot. More exactly, it takes place for the time interval before 
the moment of the Universe clearing due to hydrogen recombination. Assuming that the burst 
radiation has no effect on the medium parameters the burst size does not depend also from the 
burst power. Calculated exact profiles of intesity distribution give an opportunity to determine 
relation between angular size of a spot and physical distance to a source. Very important feature 
is also a presence of radiation polarization in the spot. Polarization is standard one for the 
case of star radiation scattered in a shell of gas with free electrons. The plane of polarization 
contains the ray of sight and direction towards the star center. Degree of polarization grows 
when moving off a spot center. However the most probable observed polarization degree will 
be about 15% on average because of fast brightness decrease to the edge of a spot. 

This research is supported by Russian Foundation for Basic Research (grant 08-02-000493). 
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Figure 1: Initial (isotropic) distribution of radiation intensity. 
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Figure 2: Distributions of the mean intensity at z — 10 (left profile) and z = 20 (right profile) 
(z = 2000). Asteriscs on the abscissa mark distances equal to conformal times of the burst for 
corresponding z. 
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Figure 3: Distributions of the mean intensity at the present epoch for different z . In the order 
of maximum intensity growth they correspond to z =3000, 2000, 1600, 1400 and 1200. In the 
same order the steepness of the leading (right) front of distribution grows. Asteriscs on the 
abscissa mark distances indicated in Tabl. 2. 
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Figure 4: The same as in the preceeding Figure but for z =1100 (1), 1200 (2) and 1400 (3). 
Asteriscs mark corresponding conformal times of the burst (see Tabl. 1). 
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Figure 5: Angular distributions of radiation intensity at z = at different distances r (see Tabl. 
1) from the burst center and for different z : thin continuous lines - z — 2000, dashed lines - 
z = 3000 and thick continuous lines - zq — f600 (for two distances - N = 4 and 2 in Tabl. 2). 
Here d is an angular distance from the burst center. 
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Figure 6: The same as in the preceeding Figure for zq = 2000 in scalar approximation (con- 
tinuous lines) in comparisson with distributions for exact description of Rayleigh scattering 
(dashed lines). 
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Figure 7: Angular distributions of radiation polarization at z = at different distances r from 
the burst center (designations are the same as in the preceeding Figure). Lowest (with the 
smallest maximum) curve corresponds to distance 13640 Mpc. 



14 



